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I.  INTRODUCTION 


Many  tyf)€S  of  artillery  projectiles  are  accelerated  in  the  barrel  by  a  sabot.  The  transfer  of  the 
accelerating  force  from  the  sabot  to  the  projectile  is  accomplished  by  a  series  of  threads  or 
grooves  in  the  later.  These  surfaces  are  referred  to  as  buttress  threads  or  grooves  in  the  follow¬ 
ing.  Aerodynamically,  these  mechanical  engagement  devices  are  roughness  elements  although 
the  roughness  heights  can  be  large  relative  to  the  boundary  layer  thickness.  They  produce  two 
effects:  one  is  an  increase  in  the  axial  force,  and  the  other  an  increase  in  the  thickness  of  the 
boundary  layer  which  can  result  in  an  interaction  with  the  rear  end  of  the  configuration.  These 
effects  may  be  significant,  hence  they  must  be  included  in  the  prediction  of  the  aerodynamic 
characteristics  of  such  projectiles.  A  shadowgraph  of  a  typical  high  speed  kinetic  energy  projec¬ 
tile  in  flight  is  shown  in  Figure  1  in  which  the  buttress  grooves  and  their  effect  on  the  boundary 
layer  are  evident. 

The  goal  of  this  work  is  to  review  prediction  methods  for  turbulent  boundary  layers  over 
rough  surfaces,  with  the  objective  of  selecting  the  ones  that  best  fit  the  ongoing  activity  of  the 
Computational  Aerodynamics  Branch,  Launch  and  Flight  Division,  US  Amiy  Ballistic  Research 
Laboratory,  (BRL)  Aberdeen  Proving  Ground,  Marydand,  and  then  to  implement  and  validate 
them. 

Most  roughness  elements  employed  in  ballistics  are  circumferential.  From  an  aerodynamic 
point  of  view,  they  are  two-dimensional.  As  a  result,  the  seiu'ch  covers,  mainly,  investigations 
concerning  two-dimensional  roughness. 

II.  CLASSIFICATION  OF  PREDICTION  METHODS 

Prediction  methods  can  be  classified  into  two  main  approache::  1)  Correlation 

methodologies;  and  2)  Analytical  predictive  techniques.  The  first  ones  are  based  on  existing 
experimental  (or  empirical)  data  and  are  useful  as  long  as  sufficient  data  are  available  under 
conditions  reasonably  close  to  the  design  problem.  It  is.  however,  of  little  help  if  the  data  are 
sparse  or  missing.  Current  activities  at  BRL  include  Mikhail's'  correlation  of  available  data  on 
rough  projectiles. 

A  classification  tree  of  predictive  methods  is  shown  in  Figure  2.  The  main  methods  for  tlie 
analysis  of  boundary  layers  are:  1 )  integral  metluxis;  and  2)  differential  methods.  Both  classes 
have  been  successfully  used  for  the  case  of  turbulent  boundary  layers  over  rough  surfaces. 
However,  the  differential  methods  can  be  modified  more  e;isily  to  handle  the  three-dimensional 
case  of  a  configuration  at  an  angle-of-attack  and  the  case  of  the  interaction  between  a  body 
boundary-layer  and  stabilizing  fins.  Also,  the  differential  methrxls  have  an  advantage  in  the 
present  work  becau.se  they  can  be  incorporated  into  existing  numerical  techniques  and  codes  in 
the  field  of  computational  aerodynamics. 


The  differential  methods  can  be  divided  into  two  major  approaches:  1)  distributed  roughness; 
and  2)  discrete  element. 

The  foundation  for  the  distributed  roughness  approach  is  the  well  established  empirical 
relationship  between  sand-grain  roughness  and  the  skin  friction  and  the  velocity  profile.  For 
roughness  other  than  that  of  sand  grain,  empirical  data  is  used  to  find  an  equivalent  sand-grain 
roughness.  This  approach  will  be  elaborated  on  in  the  next  Chapter. 

The  discrete  element  approach  estimates  the  fomi  drag  and  the  heat  transfer  to  the  individual 
roughness  elements.  These  quantities  are  added  to  the  friction  drag  and  the  heat  transfer  to  the 
smooth  .sections  between  the  excrescences.  This  approach  can  be  further  divided  into  two 
branches,  according  to  the  grid  u.sed  in  the  numerical  solution;  1)  micai  grid;  and  2)  macro 
grid. 

In  the  first  case,  the  grid  follows  the  actual  su:.*"-  -  e  Thus,  it  has  to  be  fine  enough  to  describe 
the  geometry  and  the  turbulence  quantities  in  detail.  In  principle,  this  approach  will  solve  the 
complete  (i.e.,  time  dependent)  Navier-Stokes  equations  or  the  time-averaged  ones.  This  method 
was  used  recently  by  Sahu  and  Danberg*  to  study  supersonic  flow  over  a  rotating  band.  Baysal 
and  Stallings^  studied  fiow  over  a  cavity  and  Venkataphaty,  et  al"*  considered  the  flow  over  a 
two-dimensional  groove  having  a  square  cross-.section. 

In  spite  of  the  demonstrated  success,  to  date,  this  procedure  has  had  application  only  in 
configurations  having  a  small  number  of  roughness  elements.  Many  configurations,  in  the  field 
of  ballistics,  have  large  series  of  roughness  elements.  In  these  cases,  the  computational  require¬ 
ments  are  beyond  the  present  capacity  of  codes  and  machines. 

The  second  branch,  that  of  macro  grid  ignores  the  fine  details  of  the  flow  near  the  wall.  The 
actual  surface  is  represented  by  a  smooth  one.  The  equations  of  motion  are  altered  to  account  for 
the  blocking  effect  of  the  protruberances;  the  estimated  fomi  drag  of  the  discrete  elements  is 
taken  into  account  in  the  momentum  equation;  the  estimated  Iteat  transfer  is  taken  into  account 
by  the  energy  equation.  Hodge  and  Adams^  used  this  approach  for  sand-grain  roughness  at 
supersonic  and  hypersonic  Mach  numbers.  More  recently.  Taylor,  et  al^  extended  this  approach 
for  two-dimensional,  rib-type  roughness. 

The  following  conclusions  are  made,  based  on  the  above  overview;  The  micro  grid  branch 
has  the  greatest  potentitil  tind  is  highly  recommended  because  it  requires  the  least  specific 
empirictil  data.  However,  it  is  presently  limited  to  configurations  having  only  a  smtill  number  of 
roughness  elements.  It  is  expected  that  the  application  of  th.is  brtmch  will  expand,  as  computing 
capacity  grows  in  the  future. 


The  main  issue,  for  the  analysis  of  configurations  having  many  roughness  elements  is  between 
the  distributed  roughness  approach  and  the  macro  grid  branch  of  the  discrete  element  approach. 
The  later  introduces  physical  reasoning  to  the  analysis.  Nevertheless,  it  depends  on  empirical 
data  and  contains  assumptions  that  need  verification,  especially  for  application  to  supersonic 
N;..cn  numbers  and  dense  roughness.  At  present,  it  does  not  provide  an  advantage  over  the  other 
and  more  established  methods.  Hence  the  distributed  roughness  approach  was  selected  for  the 
present  study. 

III.  SAND  ROUGHNESS 
1.  REVIEW  OF  FUNDAMENTALS 

It  is  commonly  assumed,  in  the  analysis  of  turbulent  boundary  layers,  that  Reynolds  stresses 
dominate  in  the  inner  region,  e.xcept  for  the  wall  layer.  This  assumption,  together  with  similarity 
considerations,  led  to  the  well  known  Law  of  the  Wall; 


~  =  -  /n  — 

U^  K  ‘ 


+  C  . 


Schlichting's^  values  for  the  Von  Karman  con.stant,  k.  and  the  free  constant.  C.  are  0.40  and  3.5. 
respectively.  These  values  will  be  used  in  the  present  work. 

Roughness  elements,  are  a  source  of  additional  turbulent  eddies  near  the  w-all.  These  eddies 
increa.se  the  Reynolds  stresses  in  the  wall  layer.  It  is  well  established  (e.g.  Clauser,*  Rotta*\)  that 
the  effect  of  the  additional  stresses  is  a  downward  shift  of  the  logarithmic  velocity  profile.  The 
Law  of  the  Wall  for  flow  over  rough  surface  becomes. 

—  =  -  /n—  +  C  -  —  .  (2 


The  change  of  the  intercept  depends  on  the  normalized  roughness  height,  hu.j/v,  the  geometry  of 
the  roughness  elements,  their  density  and  their  arrangement  on  the  wall. 

2.  EMPIRICAL  DATA 


The  first  and  most  (.]uoted  study  of  the  effect  of  sand  roughness  on  pipe  flow  is  that  of 
.Nikuradse.”'  He  validated  the  Law  of  the  Wail  for  rough  surfaces  and  proposed  three  regions  of 
roughness 


<  7  :  hydraulically  smooth, 

7  <  <  70  :  intermediate  rough,  (3) 

hoUr 

70  <  — ^  :  fully  rough. 

The  dependence  of  the  normalized  displacement  of  the  logarithmic  velocity  profile,  relative  to 
the  smooth  case,  is  shown  in  Figure  3.  For  the  fully  rough  region  it  was  found,  based  on 
Nikuradse’s  data,  that 


Ah  I 

=  i  /n—  -  3.0 


(4) 


Goddard,"  Fenter,'‘  and  Reda,  et  al'^  extended  the  data  base  to  supersonic  Mach  numbers. 
Their  results,  concerning  the  displacement  of  the  logarithmic  velocity  profile  are  shown  in  Figure 
3,  in  comparison  with  those  of  Nikuradse.  It  is  api)arenl  that  the  displacement  is  independent  of 
Mach  number,  provided  the  kinematic  viscosity  is  evaluated  at  the  wall. 

The  amplification  in  skin  friction,  which  is  the  ratio  of  rough  wall  to  smooth  wall  skin-friction 
coefficients,  is  shown  in  Figure  4.  The  following  relationship  has  been  found  by  the  present 
authors  to  fit  the  data. 

Cf  hjU^ 

=0.11  +  0.89  log^  .  (3) 

'“fo  '^w 

IV.  EQUIVALENT  SAND  ROUGHNESS 
1.  THE  CONCEPT  OF  EQUIVALENT  SAND  ROUGHNESS 

Three  years  after  Nikuradse"’  obtained  the  data  on  How  in  pipes  roughened  by  closely  packed 
sand  grains.  Schlichting’’^'^  published  results  of  measurements  of  the  drag  produced  by  other 
fonns  of  roughness.  His  configurations  were  made  of  spheres,  spherical  segments,  cones  and 
short  angles,  arranged  in  regular  patterns.  Some  configurations  yielded  more  drag  than  that 
produced  by  sand  roughness  of  the  same  height,  w  hile  others  gave  drag  that  w  as  smaller. 

Schlichting  proposed  to  consider  sand-grain  roughness  as  a  standard  and  introduced  the 
concept  of  equivalent  sand  roughness  which  is  the  size  of  sand  grains  that  produce  the  same  drag 


4 


as  an  actual  configuration,  under  the  same  flow  conditions.  The  range  of  the  ratio  of  the  size  of 
the  equivalent  sand  roughness  to  the  height  of  the  roughness  elements  in  his  test  varied  between 
0.15  and  3.8,  demonstrating  the  importance  of  the  shape  and  the  arrangement  in  detemiining  the 
effect  of  roughness. 


It  is  assumed  that  the  wall  velocity  profile  with  actual  roughness  is  the  same  as  that  for 
equivalent  roughness.  Introducing  the  index  ’es’  for  equivalent  sand,  the  condition  for 
Cf  =  Cf^j.  gives; 

Au  /^^u^  ,,, 


In  the  region  of  fully  rough  flow; 


i  .  D 


Hence,  the  equivalent  sand  roughness  can  be  obtained  by; 


1  ,  1  ,  *^es^'i  ^ 

-  (n—  +  D  =  -  +  D, 


Introducing  =  -3.0,  as  found  by  Nikuradse,  and  arranging: 


~  =  exp|K(3.0  D)J  .  (9) 

This  relationship  gives  a  means  of  determining  the  equivalent  sand  roughness,  given  the  change 
in  the  intercept  of  the  Law  of  the  Wall. 

2.  REVIEW  OF  EXISTING  CORRELATIONS 

The  effect  of  spacing  on  two-dimensional  roughness  uas  first  studied  by  Morris.  He 
classified  the  flow  fields  into  three  types:  skimming,  wake-inteiference,  and  isolated  roughness. 

.A  pioneering  attemr'  to  ijuantify  the  effect  of  the  spacing  of  two-dimensional  roughness 
elements  was  performed  by  Beuerinann’^'  in  l%.‘i.  Fur  roughness  composed  of  uansverse  scjuare 
bars  he  correlated  the  change  of  the  intercept  by  the  lelationship: 


D  =  17.5.5(1.634  log  Xh  ’  > 
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where  the  roughness  density  parameter  is  *he  pitch  or  wave  length  to  height  ratio,  =  p/h  . 

In  1969,  Dvorak’^  introduced  pitch  to  width  ratio  as  a  density  of  roughness  elements 
parameter,  Xg  =  p/a .  For  parameter  values  less  than  about  5.0,  and  for  square  roughness,  he 

corroborated  Bettermann’s  coirelation.  For  larger  values,  part  of  Schlichting’s  data  ana  that  of 
other  sources,  were  used  to  establish  a  comptcnviuary  correlation. 

f  17.35(1.625  log  X.,  -  1.0)  •  ^  ’ 

D  =  i  6  4/  a  jj, 

I  -5.95(1.103  log  Xg- 1.0)  .  Xg  >  4.68)  . 

Note:  The  small  difference  between  the  constant  multiplying  log  X  in  Bettetmann’s  article  and 
that  published  later  by  Dvorak,  is  of  little  practical  significance. 

The  above  relationship  between  the  displacement  of  the  velocity  profile  and  the  density  of  the 

1 8 

roughness  was  adopted  by  Cebeci  and  Smith  in  their  text  on  the  analysis  of  turbulent  boundary' 
layers. 

In  1973  Simpson'^*  suggested  a  more  general  interpretation  of  the  Bettermann  and  Dvorak 
correlation,  by  replacing  their  roughness  density  parameters  by  the  ratio  of  reference  area  to  total 
roughness  frontal  area.  The  reference  area,  S,  is  defined  as  the  smooth  surface  before  adding  the 
roughness. 


The  rational  for  this  generalization  lays  in  the  fact  that  for  fully  rough  surfaces,  the  tangential 
forces  are  predominately  due  to  form  drag  of  the  elements,  which  is  proponional  to  their  frontal 
area.  By  including  additional  data,  he  di.scovered  two  empirical  findings.  First,  for  low  values  of 
the  density  parameter  two  branches  exist,  depending  on  the  fomiation  or  nonexistence  of 
transverse  vortices  between  rows  of  the  roughness  element;;.  Second,  the  shape  of  the  roughness 
elements  is  an  important  parameter  that  should  not  be  overlooked. 

For  the  case  of  two-dimensional  roughness,  the  Simpson  interpretation  of  the  density- 
parameter  is  identical  to  that  of  Bettermann. 

The  two  branches  of  the  correlations  represent  two  types  of  flow  in  the  vicinity  of  the  rough¬ 
ness  elements.  For  small  longitudinal  spacings,  vortices  are  formed  in  the  grooves  and  the 
external  flow  skims  over  the  tops  of  the  roughness  elements  and  over  the  enclosed  vortices.  This 
kind  of  geometry  is  called  "d  type  "  because  the  vehxity  profiles,  in  pipe  flow,  scales  with  the 
diameter.  For  large  spacings.  each  element  produces  a  wake  that  interacts  with  the  following 
element.  This  is  the  "k  type  "  roughness,  where  the  height  of  the  elements  is  the  length  scale.  An 
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extreme  case  of  this  type  is  that  of  isolated  roughness  elements. 


Young  and  Patterson,^"  illustrates  the  two  types. 


Figure  5,  reproduced  from 


An  attempt  to  consider  the  shape  of  the  elements  in  a  correlation  was  reported  by  Dirling^'  in 
1973,  in  which  he  included  the  mean  windward  surface  inclination.  His  combined,  roughness 
density  and  shape  parameter  is: 


where  d  is  the  average  roughness  element  spacing  and  is  tlie  windw/ard  surface  area  wetted  by 

2'^ 

the  flow  and  A^  is  the  frontal  area  of  a  roughness  element.  In  1975,  Grabow  and  White  verified 
Dirling’s  correlation  by  showing  good  agreement  with  additional  data. 


Dalle  Donne  and  Meyer^^  experimented  with  anulli  whose  inner  tubes  were  roughened  by  ribs 
having  rectangular  cross-sections.  They  correlated  their  data  and  those  of  18  previous  research¬ 
ers,  using  the  roughness  ('ensity  parameter; 

^h  -  1,  -  ^h  h  • 


For  square  and  near  square  rods  (0.95  <  h/a  ^  1.05),  their  correlation  gives: 


D  = 


,>•<1.73 

5.5  -  9.3  Xj, 


+114(1 

5.5  -  1.04  X|, 


1.0  ^  <  6.3  , 

6.3  <  Xi*  i  160.0  . 


(12) 


TABLE  1  contains  a  list  of  the  roughness  density  paranieieis  used  in  the  past.  Note  that  for 
square  bars,  nomial  to  the  mainstream,  all  the  lust  four  p;iiametcis  are  identical  and  the  Dalle 
Donne  and  Meyer  ptirametei  is  related  to  them  by  X|,  =  Xj^  1.0  . 
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TABLE  1.  Definitions  of  Roughness  Density  Parameters 


Authors 

Year 

Parameter 

Reference 

Bettemiann 

1965 

^h 

n 

16 

Dvorak 

1969 

^a 

=  £ 
a 

17 

Simpson 

1973 

S_ 

Sf 

19 

Dirling 

1973 

A 

- 

^  J 

21 

Dalle  Donne  &  Meyer 

1976 

^h 

P  •  a  ,  a 

'  h  "  •  i; 

23 

•^s 


If) 


Present 


1988 


3.  A  NEW  CORRELATION 


Since  the  above  mentioned  correlations  have  been  published,  additional  data  have  been 
reported.  Furthermore,  Coleman,  Hodge  and  Taylor^'*  re-evaluated  Schlichting’s  experimental 
data.  Their  corrected  values  for  the  equivalent  sand  roughness  are  smaller,  and  in  some  cases 
much  smaller,  than  those  reported  in  Reference  14  and  used  in  the  previous  correlations.  Refer¬ 
ence  24  re-evaluated  three  out  of  four  of  the  groups  of  roughnesses  tested  by  Schlichting.  The 
corrected  equivalent  sand  roughness  for  the  fourth  group,  that  of  short  angles,  was  obtained  using 
a  correlation  of  the  available  corrected  versus  uncorrected  values.  The  uncertainty  involved  in 
this  approximation  is  smaller  than  the  scatter  in  the  data. 

The  study  of  the  available  data,  including  the  corrected  data,  showed  that  only  Simpson's 
roughness  density  parameter  represents  the  data  for  three-dimensional  roughness  configurations 
close  to  that  of  two-dimensional  ones.  Also,  the  inclusion  of  a  shape  factor,  as  proposed  by 
Dirling,  is  needed  in  order  to  correlate  data  from  roughness  elements  of  different  shapes.  The 
new  roughness  density  parameter  is; 


A 


s 


if  ^ 

^f  Y'^S 


1  /> 


and  the  correlation  is  shown  in  Figure  6.  The  sources  of  the  data  are  summarized  in  TABLE  2. 
Note  that  the  previous  correlations,  plotted  in  Figure  6,  arc  only  valid  for  square  or  near  .square 
bars.  In  all  these  cases,  the  shape  factor  becomes  unity  so  that  the  former  roughness  density 
parameters  are  all  identical  to  the  current  one.  The  power  of  the  shape  factor  in  the  current 
parameter  has  been  adjusted  so  that  it  correlates  Schlichting’s  corrected  data  for  all  four  groups 
of  roughness  elements.  In  the  case  of  circular  rods  (References  32  and  33),  the  streamwise  area, 
below  the  point  of  maximum  width,  was  accounted  equal  to  the  frontal  area 
(A,/Ap  =  0.5 -t- n/2.0). 


The  correlation  of  the  three-dimensional  data  is  consistently  lower  than  that  of  two- 
dimensional  roughnesses.  This  difference  is  expected  because  the  end  effects  of  the  finite  width 
elements  reduce  their  fomi  drag  relative  to  continuous  transverse  ones. 

The  model  is  biased  toward  data  obtained  on  plane  trough)  surfaces.  For  low  values  of  the 
roughness  density  parameter,  the  Bettennann  relationship  of  Equation  (10),  which  is  based  on 
parametric  wind  tunnel  tests,  is  retained.  However,  the  large  spread  in  the  later  data  adds  uncer- 
iainty,  w'hich  only  can  be  solved  by  the  addition  of  systematic  new  data.  For  moderate  and  liigh 
values  of  the  density  parameter,  this  work  proposes  for  two-dimensional  roughness: 


TABLE  2.  Sources  of  Data  Used  in  the  Correlation 


Authors 

Reference 

range 

Comments 

Schlichting 

14 

0.9  -  9.7 

spheres 

2.6-  15.3 

spherical  segments 

5.3-  10.7 

cones 

6.7  -  13.3 

short  angles 

Bettemiann 

16 

2.65  -  4.18 

low  speed  wind  tunnel, 
square  bars 

Dalle  Donne  &  Meyer 

23 

4.08  -  61.5 

annular  tube,  rough  inner 
rod.  Correlation  includes 

data  of  several  researchers 
and  cover  2.0  <  3.),  <  16(,).() 

Webb,  et  al 

25 

10.0  -  40.0 

tubes 

Han,  et  al 

26 

5.0-  15.0 

parallel  plates,  only  data  for 
transverse  ribs  used 

Liu,  et  al 

27 

2.0  -  96.0 

water  tunnel 

Pincau,  et  al 

28 

4.0  ^ 

Perry  &  Joubert 

29 

4.0 

low  speed  wind  tunnel, 

)  square  bars 

Antonia  &  Luxton 

30 

4.0 

Antonia  &  Wood 

31 

2.0  y 

Furuya,  et  al 

32 

2.0-64.0 

transverse  circular  rods 

Sherif  &  Giimley 

33 

10.0 

10 


(13) 


D  =  < 


17.35(1.634  log  As  -  1.0) 

2.2 

9.55(1.0  -  0.686  log  As) 


1.4  <  As  <  4.89  . 

4.89  <  As  <  13.25  , 
13.25  ^  As  s  100.0. 


For  dense  arrangements  of  two-dimensional  roughness  elements,  vortices  form  in  the  cavities 
and  the  external  flow  skims  over  the  tops  of  the  elements  and  tlie  vortices.  In  this  case,  the  shape 
factor  may  be  different  from  that  for  sparse  roughness  elements.  Since  no  additional  pertinent 
data  were  found  for  this  region,  the  left  branch  of  the  correlation  (Equation  ( 13))  should  only  be 
used  for  square,  or  near  square,  bars.  The  scarcity  of  data  for  tlirec  dimensional  roughness  at  low 
and  intennediate  values  of  the  density  parameter  prevents  a  general  mcxlel  for  this  case.  A  partial 
model  for  three-dimensional  roughness  is: 

D  =  9.7(l.(l  -  0.794  log  As)  ,  16.0  <  A^  <  200.0  .  (14) 

4.  THE  EQUIVALENT  SAND  ROUGHNESS 

Using  Equation  (9).  the  equivalent  s.uhI  roughness  for  two-dimensional  roughnesses  becomes: 


^es 

h 


().()()3215  As-"’’- 

^  8.0 


y  1.51.71 


1.4  As  ^  4.89 
4.89  <  As  <  13.25 
13.25  S  As  <  100.0 


(15) 


for  three-dimensional  roughnesses: 


160.77  Aj, ' 


16,0  <  A^  <  20().()  . 


(16) 


These  correlation  equations  are  presented  in  Figure  7.  The  present  left  branch  is  identical  to  that 
of  Bettermaim  and  very  close  to  that  of  Dirling.  On  the  oilier  liaiul.  the  present  right  brancii  is 
much  higher  than  the  previous  correhitions  by  Dvorak  and  Dirling.  Due  to  the  variance  in  l' 
density  parameters,  comparison  of  the  equivalent  sand  roughness  is  only  possible  for  the  case 
square  bars. 
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5.  DEPENDENCE  ON  MACH  NUMBER 


As  discussed  in  Chapter  HI,  the  displacement  of  the  logarithmic  velocity  and  the  gain  in  skin 
friction  for  sand  roughness  is  independent  of  Mach  number,  provided  the  Reynolds  number  of 
the  grain  size  is  evaluated  at  wall  conditions.  Voisiner"*  reached  the  same  conclusion  for 
screens. 

On  the  other  hand,  the  only  datum  found  for  two-dimensional  roughness  at  a  supersonic- 
speed,  is  that  of  Berg,^^  which  shows  an  effective  sand  roughness  that  is  lower  than  predicted  by 
the  low  speed  correlations.  The  results  concerning  skin  friction  and  the  intercept  of  the  logarith¬ 
mic  velocity  are  consistent,  indicating  high  cjuality  data.  Nevertheless,  since  this  is  a  single 
datum,  no  reliable  relationship  between  the  equivalent  sand  roughness  and  Mach  number  can  be 
devised.  Such  a  relationship  will  have  to  be  based  on  calibration  of  computations  of  the  drag  of 
roughness  configurations  for  which  data  are  available. 


V.  APPLICATION  IN  COMPONENT  BUILD-UP  METHODOLOGIES 
1.  COMPONENT  BUILD-UP  IN  DRAG 

The  new  correlation  can  be  used  for  an  engineering  type  estimate  of  the  increment  in  friction 
drag  due  to  surface  roughness.  The  results  of  such  estimates  can  be  used  with  Component 
Build-Up  (CBU)  methodologies,  because  they  add  on  the  various  contributions  to  the  drag. 

A  schematic  of  the  analysis  is  shown  in  Figure  X.  Using  the  new  correlation,  the  equivalent 
sand  roughness,  ^or  a  given  roughness  geometry  is  evaluated.  Also,  for  given  flight  conditions 
and  size  of  the  flight  vehicle,  the  smooth  wall  average  skin-friction  coefficient  is  estimated. 
Knowing  these  two  quantities,  the  gain  in  skin  friction,  Cf/Cf^  ,  due  to  roughness,  is  evaluated. 
The  additional  axial-force  coefficient  is 


^^A  “  Sr  ’ 


(17) 


where  ACf  =  Cj-  -  Cf^  . 

2.  LONG-ROD  PROJEtTflLE  TEST  CONFIGUR.ATIONS 


The  component  build-up  scheme  and  the  Navier-Stokes  numerical  computation,  to  be  consid¬ 
ered  in  subsequent  Chapters,  were  checked  by  applying  them  to  the  wind-tunnel  models  tested  by 
Brandon  and  Von  Wahlde.'^^  The  models  were  tested  in  the  Naval  Surface  Weapons  Center 


supersonic  wind  tunnel  and  consisted  of  a  Sears-Haack  nose,  a  forebody,  an  afterbody  and  a  fin 
section  as  shown  in  Figure  9.  Several  forebixlies  and  afterbodies  were  constructed  in  different 
lengths  so  that  the  overall  length  could  be  varied  between  20  and  35  calibers.  The  forebody 
sections  were  either  smooth  or  grooved  as  shown  in  Figure  10.  The  afterbody  sections  were 
either  smooth,  grooved  or  covered  with  standard  machine  threads  {30  threads/caliber).  Thus  a 
number  of  different  kinds  of  roughness  characteristics  could  be  produced  to  simulate  typical 
saboted  projectiles.  Tests  were  performed  with  and  without  the  fins.  The  wind  tunnel  model  has 
an  abrupt  6.2%  decrease  in  diameter  at  the  fore-  to  aft-body  junction. 

The  wind  tunnel  tests  of  the  models  were  ctirried  out  over  a  Mach  number  range  of  3.5  to  5.0. 
However,  the  majority  of  the  tests  were  focused  on  Mach  number  5.0.  The  Reynolds  number  in 
the  wind  tunnel  was  relatively  small  compared  to  free-flight  conditions,  particularly  at  the  high 
Mach  numbers. 

3.  COMPONENT  BUILD-UP  RESULTS 

The  roughness  density  parameter  for  these  buttress  grooves  is  small  enough  to  be  on  the  left 
branch  of  the  new  correlation.  According  to  a  comment  in  Section  IV-3.  the  shape  factor  should 
not  be  used  in  this  case.  Thus,  the  equivalent  sand  roughness  was  estimated  as  if  the  cross- 
section  of  a  groove  was  a  square.  The  ratio  h^Th  was  found  for  two  cases.  The  first,  for 

=  p/h  and  the  other  for  account  the  width  of  the  cavities  at 

their  top  end.  An  average  value  of  the  two  cases  is  used.  For  the  (plain)  threads,  the  data  of 
Fenter  for  v  grooves  was  used.  The  estimated  values  are: 


fomi 


h,/h 


Buttress  Grooves  0.52 

Threads  0.55 


The  smooth  wall  skin  coefficient  was  obtained  from  the  NSWC  Aeropredictioiv^^  code.  For 
configurations  with  mixed  roughness  {buttress  threads  on  the  forebody,  threads  on  the  afterbKxly), 
separate  estimates  were  perfomied  for  each  section.  The  analysis  followed  the  scheme  described 
above  and  the  results  arc  show'ii  in  Figure  I  1  in  comparison  with  wind-tuntiel  data,  bttsed  on 
Appendix  D  of  Referetice  36.  In  more  than  half  the  cases,  the  agreement  between  prediction  and 
test  data  is  gocxl.  In  the  other  cases,  the  deviations  do  not  show  a  consistent  trend  and  are  of  the 
same  order  of  magnitude  as  the  scatter  in  the  data. 
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VI.  ALGEBRAIC  TURBULENCE  MODELS  FOR  NUMERICAL  CODES 


In  Sections  IV  and  V,  the  conccpv  of  equivalent  sand  roughness  was  explored  in  some  detail. 
The  added  drag  was  estimated  for  a  given  configuration  using  the  component  build-up  technique. 
Another  approach,  that  is  potentially  more  general,  is  to  modify  the  algebraic  turbulence  models 
used  in  existing  numerical  Navier-Stokes  solvers  to  account  for  surface  roughness.  The  distrib¬ 
uted  roughness  approach  does  not  consider  the  details  of  the  roughness  elements,  so  it  is  not 
simply  a  matter  of  changing  local  boundary  conditions  at  the  surface.  The  effect  of  the  rough¬ 
ness  is  incorporated  by  changing  the  equations  which  define  the  relationship  between  the  tur¬ 
bulent  eddy  viscosity  and  the  local  flow  variables.  Since  most  codes  use  some  fomi  of  a  two- 
layer  mixing-length  method  for  computing  the  eddy  viscosity,  the  effect  of  roughness  is  to 
modify  the  mixing-length  fomiulation.  Several  researchers  have  studied  this  problem  and  the 
results  of  their  analyses  are  discussed  in  this  section  followed  by  details  of  the  incorporation  of 
one  of  them  into  a  parabolized  Navier-Stokes  code  and  its  application  to  the  Brardon  experiment. 

1.  VAN  DRIEST  MODEL 

The  first  mixing-length  model  that  considers  wall  roughness  was  developed  by  Van  Driest^*^ 
(the  superscript  +  indicates  wall  units): 


/*  =  Ky'^ 


exp 


^  -'P  2-g; 


,  h;  <  60 


(18) 


The  value  h^  =  60  marks  the  beginning  of  the  region  of  fully  rough  flow.  For  this  value  the 
wail  damping  is  eliminated  and  the  mixing  length  becomes: 


r  =  tcv^  at  hr  =  60 


According  to  this  model,  /  vanishes  at  the  wall  in  the  entire  range  i  f  applicability,  as  illustrated  in 
Figure  12.  The  initial  slopes  of  the  T  vs  y*  curves  are  equal  to  <,  except  for  the  case  of  a 
smooth  surface,  where  the  initial  slope  vanishes.  To  apply  this  method  to  larger  roughness 
(hej;  >60)  requires  additional  empirical  information  concerning  the  displacement  in  y  of  the 
effective  wall. 


2.  HEALZER  MODEL 

Several  researchers  proposed  extensions  of  the  mixing-length  model  into  the  region  of  fully 
rough  boundary  layer,  by  allowing  the  mixing  length  to  exceed  Ky*  and  become  non  zero  at  the 
wall.  The  one  by  Healzer,^*^  et  al  was  considered  for  application  because  it  w-as  successfully  used 
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by  Lin  and  Bywaier.'*^  It  has  been  selected  for  inclusion  in  the  PNS  code  because  it  is  compat¬ 
ible  with  the  Baldwin-Lomax  turbulence  model  already  in  the  code,  it  is  directly  adopted  to 
compressible  flow,  and  it  is  not  limited  in  roughness  height  which  can  be  considered. 


For  intermediate  roughness,  the  mixing  length  and  it;  inl*ial  slope  vanish  at  the  wall.  For  a 
fully  rough  wall,  t*  equals  at  the  wall.  The  wall  value  increases,  as  h*  increases,  as  shown  in 
Figure  1 3. 

3.  HAN  MODEL 

More  recently,  Han"*'  developed  a  new  formulation  for  the  mixing  length  for  smooth  surfaces 
and  extended  it  to  rough  ones.  His  approach  is  different  from  all  previous  ones  in  relating  the 
mixing  length  to  the  velocity,  rather  than  to  the  distance  from  the  wall.  Four  models  were 
studied  by  Han,  the  preferred  one  being 

r  ~  R^[exp(KU*)  -  exp(-KU''))  ,  (20) 

with  E  =  9.025  and  R.  the  roughness  amplification  factor,  a  function  of  the  shape  of  the  rough¬ 
ness  and  of  h’’.  For  sand  grain  roughness: 

R  =  0.3036  h;  .  h;  >  100  .  (21) 

For  h*  <  1(K),  Han  provides  tabulated  values  of  R.  Study  of  the  data  showed  that  the  above 
expression  can  be  used  down  to  h,  =  50  with  a  very  small  inaccuracy.  The  following  model  is 
proposed 

r  0.3036  h;  ,  h;  >  50  . 

R  =  (22) 

i  0..M)36h;  -  0.043(50  -  h")  ,  20  <  h!  <  50  . 
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It  should  be  noted  that  the  Han  model  is  based  on  the  incompressible  logarithmic  velocity  profile 
and  as  a  consequence  it  is  not  apparent  how  this  formulation  can  be  extended  to  high  speed 
compressible  flows. 


4.  THE  OUTER  LAYER 


Clauser*^  reached  a  conclusion  that  the  turbulence  level  and  the  Reynolds  stresses  in  the  outer 
layer  are  not  affected  by  wall  roughness.  Since  then,  most  researchers  have  accepted  this  find¬ 
ing.  Thus  in  the  following,  the  Baldwin-Lomax  outer-flow  computational  procedure  is  used. 


Vll.  APPLICATION  TO  PARABOLIZED  NAVIER-STOKES  CODE 
1.  PARABOLIZED  NAVIER-STOKES  CODE 


The  Paralxjlized  Navier-Stokes  equations  (PNS)  have  been  solved  using  a  technique  devel¬ 
oped  by  Schiff  and  Steger.*^^  These  equations,  in  strong  conservation  form,  and  simplified  using 
the  thin  shear  layer  approximation  can  be  written  in  generalized  coordinates  as  follows: 


A  A  A 

"Es  .OF  50  ^  1  3S 


(23) 


where  4.  fi.  4  generalized  coordinate  variables. 

4  =  4(^)'  's  the  longitudinal  (maa'hing)  coordinate 
n  =  ri(x,y,z),  circumferential  coordinate 
4  =  near  nomial  coordinate. 

AAA  A 

The  inviscid  flux  vectors  Ej;,  F,  and  G  and  the  viscous  matrix  terms,  S,  are  functions  of  the 
dependent  variables,  q; 


A 

M 


=  J  ' 


pu 

pv 


(24) 


where  p  is  the  density,  pu,  pv.  and  pw  are  the  mass  fluxes  in  the  three  coordinate  directions  and  e 
is  the  total  energy  per  unit  volume  The  local  pressure  is  deiemiineci  by  using  the  relation: 
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P  =  (Y  -  l)le  -  0.5p(u‘  +  v‘  +  w‘)] 


(25) 


where  Y.  is  the  ratio  of  specific  heats.  The  parabcliration  of  the  equations  is  accomplished  in  the 

A 

Ej  flux  vector  by  ensuring  that  the  static  pressure  is  maintained  constant  across  the  subsonic- 
layer  near  the  wall.  Thus,  in  supersonic  free-stream  flows,  these  equations  can  be  solved  by 
marching  in  the  longitudinal  direction. 

An  initial  plane  of  data  is  required  to  begin  the  marching  and  this  is  obtained  from  the  conical 
starting  procedure.  The  starting  procedure  consists  of  marching  one  step  downstream  and 
updating  the  previous  step  by  e.xtrapolating  back  using  a  conical  flow  assumption.  This  proce¬ 
dure  is  repeated  until  a  converged  solution  is  obtained  in  which  the  chatige  in  density  between 
successive  iterations  is  less  than  10  '*'  times  the  free-stream  value. 

The  numerical  algorithm  is  formulated  using  an  approximately  factored,  implicit,  finite 
difference  scheme  developed  by  Beam  and  Warming.'*^  Fitting  of  the  outer  bow  shock  wave  has 
been  perfomied  in  these  calculations  consistent  with  the  implicit  boundary  method  fomiulated  by 
Rai  and  Chaussee."^ 

2.  IMPLEMENTATION  OF  THE  TURBULENCE  MODEL 

The  PNS  code  employs  the  Baldwin-Lomax"’*'  smcK>th  wall  turbulence  model  and  the  viscous 
flow  is  assumed  turbulent  from  the  sttuiing  solution.  The  Baldwin-Lomax  model  is  a  two  layer 
model  in  which  the  inner  layer  eddy  viscosity  is  defined  by: 

Pj  =  p/‘lci)l  (26) 


where 


/  =  ky 


and  |(i)l  =  magnitude  of  the  vorticity 


and  for  attached  flow,  the  vi.scosity  in  the  outer  layer  is  given  by: 


Pt  ~  *^^cpP^'m^*^'m^^ 


where  F(y)  =  yieui  1  -  e.\p  - 


[' 


A" 


(27) 


and  F(ypp,)  =  maximum  of  F(y)  and  Vp,  is  the  normal  distance  to  the  maximum  in  F(y). 
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K  =  0.0168 


1.6 


(Klebanoff  imemiittency  factor) 


^kleb 


0..V 


The  inner  viscosity  is  used  at  each  step  nomia!  to  the  wall  until  it  exceeds  the  outer  viscosity 
where  upon  the  outer  value  is  used. 

The  original  coding  cif  the  Baldwin-Lomax  technique  evaluated  the  eddy  viscosity  between 
grid  points  and  interpolated  to  form  the  viscosity  at  the  grid  points  for  the  subsequent  calcula¬ 
tions.  This  unnecessary  step  has  been  eliminated  by  Weinacht'^^  and  it  is  his  version  that  has 
been  used  here. 


The  Healzer,  et  al.^'^  model,  as  well  as  the  cMher  mcxiels  discussed  in  Section  VI.  is  assumed  to 
only  affect  the  inner  layer  mixing  length  and  that  the  outer  layer  formtiMtion  remains  unchanged. 
When  the  roughness  model  was  implemented,  the  solution  at  the  rou;;hness  oscillated  strongly 
and  in  many  cases  diverged.  It  was  observed  that  in  a  distributed  roughness  method,  the  rough¬ 
ness  is  viewed  as  a  change  in  the  wall  boundary  conditions  whereas  the  Healzer  model  not  only 
affects  the  wall  conditions  but  actually  changes  the  eddy  viscosity  formulation  across  the  inner 
boundary  layer.  Physically,  diffusion  of  wall  conditions  should  be  limited  to  a  zone  of  influence, 
at  least  roughly  corresponding  to  the  propagation  of  a  Mach  wave  across  the  layer.  A  convenient 
implementation  t  f  this  idea  involves  freezing  the  eddy  viscosity  at  its  value  at  the  stan  of  the 
roughness  zone  above  a  Mach  wave  originating  at  the  start  of  the  roughness.  Below  the  Mach 
wave,  the  eddy  viscosity  is  calculated  by  the  Healzer  model  but  multiplied  by  a  transition  factor 
which  smoothly  connects  the  roughness  eddy  viscosity  with  the  viscosity  above  the  wave.  A 
further  simplification  is  made  by  calculating  the  Mach  line  liKation  ba.sed  on  the  free-stream 
Mach  wave  direction  and  ignoring  the  actual  Mach  number  distribution  including  the  subsonic 
layer.  Thus  the  eddy  viscosity  distribution  is  calculated  over  the  roughness  .section  by: 

Pj  =  lij  +  (p|  -  p'j)  G  (28) 


w'here:  Pj  =  p.'ly) 

=  the  smooth  wall  eddy  viscosity  at  the  start  of  the  roughness 


p[  =  rough  wall  eddy  viscosity 


IS 


G  =  transition  factor 


The  factor,  G,  is  chosen  to  have  the  following  form. 


(  1 

=  (X  -  X(.j)  tan 

arcsinj 

i  VI 

_ 

=  location  of  Macli  wave 
Xpj  -  longitudinal  start  of  the  roughness  area 
y  =  distance  nomial  from  the  wall 

y,^,  =  Baldwin-Lomax  normal  distttnee  to  maximum  in  moment  of  voilicity. 


Note  that  when  is  significantly  larger  than  the  viscous  layer,  the  eddy  viscosity  is  computed 
entirely  based  on  the  rough  wall  model.  The  form  and  extent  of  the  transition  region  is  quite 
arbitrary  except  that  eddy  viscosity  .should  smoothly  change  from  the  smooth  wall  to  rough  wall 
value  over  as  small  of  a  region  of  the  rough  surface  as  is  needed  tor  numerical  stability. 

The  downstream  transition  from  the  rough  wall  to  the  smooili  wall  potentially  pieseiits  a 
similar  problem  of  the  discontinuous  change  In  the  ecpiaiinns  used  to  compute  the  mixing  lengili. 
A  transition  may  be  expected  where  the  smooth  wall  turbulence  diffii.;es  out  from  the  cliange  in 
the  wall  boundary  condition.  In  the  present  case,  it  was  not  necessary  to  introduce  a  transition 
process,  altnciugh  the  skin  fiiciion  decreased  to  very  low  values  behind  the  nnigliness  and  could 
possibly  go  to  zero  or  becc'me  negative  fc>r  larger  roughnesses.  Tlie  PNS  code  would  nni  be  valid 
under  these  conditions,  nor  would  the  Ilealzer  turbulence  model  be  applicable. 


3.  APPLIC/^TION  TO  A  LONG-ROD  PROJECTILE 


As  previously  noted,  the  computations  with  the  loughness  model  have  been  applied  to  the 
flow  conditions  and  configuration  tested  by  Brandon  and  Von  Wahlde.  The  model  configura¬ 
tion  is  shown  in  Figures  9  and  10. 

The  blunt  tip  of  the  Seais-Haack  nose  (0.035  caliber)  is  replaced  by  a  14.972  degree  conical 
lip  which  increased  the  total  length  to  20.757  calibers  for  the  smallest  L/D  model.  No  forces  are 
calculated  in  the  conical  tip  and  thus  the  axial  forces  are  somewhat  less  than  would  be  obtained 
in  the  experiment.  The  origin  of  the  computational  coordinate  system  is  at  the  vinual  origin  of 
the  conically  extended  nose.  The  change  in  diameter  between  the  fore-  and  aft-body  is  imple¬ 
mented  in  the  computations  as  a  conical  section,  of  about  one  caliber,  in  order  to  reduce  numeri¬ 
cal  oscillations  due  to  the  change  in  geometry.  The  PNS  computations  were  made  corresponding 
to  a  Mach  number  of  5.0,  Reynolds  number  of  18.6  x  10''  per  meter  and  a  nominal  wall  tempera¬ 
ture  of  300  K. 

The  geometry’  of  the  groove  roughness  elements  are  shown  in  Figure  10.  The  test  configura¬ 
tion  computed  here  consists  of  groove  roughness  elements  extending  over  the  entire  fore-  and 
aft-body  and  just  on  the  fore-body.  The  groove  depth  is  0.05.^5  ctilibers  and  its  pitch  is  0. 1.^34 
calibers.  The  front  face  of  each  element  is  inclined  at  approximately  a  45  degree  tingle.  As  has 
been  previously  pointed  out,  this  roughness  is  such  that  the  roughness  density  correlation  is  not 
valid  and  the  equivalent  sai'.d  roughness  height  has  been  estimated  in  Section  V  to  be  0.0278 
calibers.  However,  the  computations  have  been  performed  over  a  range  of  values  of  equivalent 
roughness  height. 

4.  PNS  RESULTS 

Figure  14  shows  the  skin  friction  coefficient  distribution  for  the  case  of  an  equivalent  rough¬ 
ness  height  of  0.0278  calibers.  The  solid  line  is  the  smooth  body  ca.se  and  the  dashed  line  is  w'ith 
groove  roughness  over  the  fore-  and  aft-body.  There  is  considerable  oscillation  at  each  change  in 
the  body  shape  which  is  frequently  observed  in  PNS  SKin  friction  results.  Skin  friction  is  an 
extre'"^ly  sensitive  result  because  it  mainly  depends  on  the  derivative  of  the  velocity  piofile. 
Increased  numerical  damping  could  possibly  reduce  this  to  some  extent  but  potentially  at  .some 
loss  in  acci'-acy. 

.•\t  the  leading  edge  of  the  roughness  there  is  a  very  sharp  peak  in  the  wall  friction  but  the 
effect  of  this  on  the  properties  of  the  boundary  layer  is  limited  by  the  outward  diffusion  as 
previously  discussed.  Downsiieam  of  the  transient,  the  friction  coefficient  becomes  netirly  a 
constant  factor  of  1 .5  higher  than  the  smooth  wall  case.  The  abrupt  transition  from  rough  to 
smooth  downstream  of  the  aft  body  induces  strong  oscillations  in  the  friction  but  these  damp  out 
after  about  four  calibers.  The  friction  considerably  under-shoots  the  return  to  the  smooth  ca.se 


and  slowly  approaches  it  from  below.  The  conical  transition  section  between  the  two  different 
diameter  cylinders  shows  some  relatively  strong  effects  of  the  changing  local  flow  on  the  rough¬ 
ness  computation. 

Figure  15  shows  the  effect  of  the  roughness  on  the  surface  pressure  distribution  is  very  weak 
under  these  conditions.  There  is  a  .slight  increase  in  pressure  near  the  start  of  the  roughness  and  a 
slight  decrease  downstream  of  the  roughness.  These  changes  reflect  the  increased  rate  of  growth 
of  the  boundary  layer  thicknesses  over  the  rough  portion  of  the  body.  Note  that  the  oscillations 
in  the  solutions  arc  much  weaker  in  the  pressure  as  compared  to  those  in  the  friction. 

Figure  16  illustrates  the  effect  on  the  velocity  profiles.  The  most  distorted  profile  occurs  at 
the  most  downstream  point.  Except  for  that  one  profile,  the  rest  show  the  expected  growth  in  the 
boundary  layer  thicknc.ss. 

The  total  forebody  drag  results  are  shown,  in  Figure  17,  as  a  function  of  equivalent  roughness 
height.  Two  cases  are  shown,  the  first  has  roughness  over  both  the  fore-  and  aft-body  cylinders. 
The  second  case  is  for  the  situation  with  roughness  only  on  the  fore-body.  The  smooth  body 
calculation  is  in  good  agreement  with  the  experiment  (it  is  four  percent  lower  because  of  the 
omitted  no.se  blunting).  There  is  a  break  in  both  curves  which  occurs  at  the  point  where  the 
Healzer  nuxlel  changes  from  h*  less  than,  to  greater  than  55.0.  The  horizontal  dotted  lines  show 
the  measured  drag  for  the  two  cases.  The  equivalent  sand  roughness  is  0.075  inches  for  the  case 
where  both  surfaces  are  roughened  and  0.1)5  inches  for  the  rough  fore-body.  Both  of  these  values 
are  significantly  different  than  the  predicted  equivalent  roughness  height.  However,  the  ex¬ 
perimental  uncertainty  in  the  drag  due  to  roughness  is  potentially  quite  large  and  this  could  have 
a  significant  effect  on  the  equivalent  sand  roughness  height  detemiined  from  Figure  17. 

XI.  CONCLUSIONS 

Based  on  a  review  of  previous  correlations  of  the  effect  of  roughness  on  the  turbulent  bound  ¬ 
ary  layer,  a  new  roughness  density  parameter  is  proposed.  It  is  a  product  of  reference  surface  to 
total  frontal  area  ratio  and  a  shape  factor  which  takes  into  account  the  average  inclination  of  the 
windward  surface  of  the  roughness  elements.  The  new  parameter  is  used  to  correlate  the  dis¬ 
placement  of  the  logarithmic  wall  profile  for  two-dimensional  and  three-dimensional  roughness 
elements.  The  equivalent  stind  roughness,  based  on  the  new  correlation,  is  identical  to,  or  in 
gotxl  agreement  with,  previous  ones,  for  small  values  of  the  roughness  parameter.  However,  for 
large  values  of  this  parameter,  the  differences  are  considerable. 

For  dense  two-dimensional  roughness,  die  ratio  of  the  size  of  the  equivalent  roughness  to 
height  of  the  actual  element,  strongly  depends  on  the  roughness  density  parameter.  However, 
because  of  the  variability  of  the  limited  avtiilahle  d;ita,  the  correiaiion  results  for  values  of  the 
roughness  parameter  less  than  five  are  restricted  to  rectangular  elements. 
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Additional  data  are  needed  in  order  to  extend  and  further  substantiate  the  correlation.  For  the 
two-dinicnsional  case,  data  for  chamfered  elements  arc  desired,  and  for  three-dimensional 
roughness,  elements  of  various  aspects  ratios  (spanwise  dimension  to  height  ratio)  in  different 
densities  and  arrangements.  The  additional  data  are  required,  mainly,  in  order  to  extend  the 
correlation  to  small  and  intermediate  values  of  the  roughness  density  parameter.  In  addition, 
compressibility  effects  at  supersonic  Mach  numbers  are  .still  poorly  understood, 

A  modified  form  of  Bettermann's  correlation  has  been  used  in  an  engineering  type,  compo¬ 
nent  build-up,  analysis  for  the  estimation  of  the  additional  axial-force  due  to  grooves  and  threads. 

A  turbulence  model,  based  on  that  of  Healzer.  et.  al,  and  the  equivalent  sand  roughness 
correlation  has  been  incorporated  in  the  Schiff  and  Steger  Parabolized  Navier-Stokes  code.  Tliis 
code  hus  been  applied  to  the  long-rod,  wind  tunnel  model  of  Brandon  and  Von  Wahlde.  It  is 
observed  that  the  Healzer  model  has  to  be  modified  to  include  a  transition  zone  when  going  from 
a  rough  to  a  smootli  surface  and  vice  verse.  This  is  required  to  obtain  stable  numerical  solutions 
and  to  account  for  the  propagation  of  information  across  the  viscous  layer  at  high  speeds.  The 
propagation  is  approximated  by  restricting  the  effect  of  the  changing  wall  conditions  to  a  region 
below  a  Mach  wave  originating  at  the  change  in  wall  conditions.  The  wave  has  been  estimated 
using  the  free-stream  Mach  number. 

The  modified  PNS  cixlc  has  been  applied  to  the  wind  tunnel  tests  of  Braiulon  and  Von 
Whalde.  When  the  predicted  equivalent  roughness  height  of  0.()27H  calibers  was  used  in  the 
code,  the  skin  friction  coefricieiii  is  increased  on  the  roughened  surface  over  the  snuxnh 
wall  case.  This  is  consistent  with  the  component  “itiild  up  method.  Additionally,  th.e  PNS  results 
pred  .  the  skin  friction  to  he  below  (up  to  20%)  the  snuxith  wall  value  on  the  surfaces 
downstream  of  the  roughness.  The  component  build-up  technique  did  not  account  for  this 
reduction. 

The  total  drag  computed  from  the  modified  PNS  code  using  the  predicted  equivalent  rough¬ 
ness  height  is  significantly  lower  than  the  drag  measured  in  the  wind  tunnel.  A  larger  equivalent 
roughness  height  corresponding  to  the  actual  rougluiess  gave  bettci'  results. 
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PREDICTION  METHODS 


Figure  2.  Classification  of  predictive  methtxis  ftx  turhulenl  lx>undary-!ayer  over  rough  surtaces  . 
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rouehness  on  the  Law  of  the  Wall  intercept 
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Figure 


k'  TYPE  ROUGHNESS 


’d*  TYPE  ROUGHNESS 


Figure  5.  Illustration  of  ’k’  and 'd'  roughnesses  (after  Reference  2U) 
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Correlation  of  the  effect  of  roughness  density  on  equivalent 
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gure  8.  ^ 
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f-itiure  14.  f’NS  c<xie  preiliclioii  r 
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Figure  15.  PNS  ccxie  preJiction  of  the  pressure  distribution  For  the  Brandon,  et  al  wind  tunnel  model. 


ta(ive  velociiv  profiles  over  the  rough  surface. 
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Ar 

As 


C 
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c  c 
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Cf 

Cfo 

d 

D 
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A  A 

..F,G 


G 

h 


width,  in  mainstream  direction,  of  roughness  element 
frontal  ai-ea  of  a  roughness  element 

parameter  in  Healzer  model 
windward  area  of  a  roughness  element 
free-stream  velocity  of  sound 
a  constant  in  the  Law  of  the  Wall 
axial-force  coefficient 

constraints  in  Baldwin-Lomax  turbulence  model 
skin  friction  coefficient 
smooth  wall  skin  friction  coefficient 
average  spacing  of  roughness  elements 

a  term  in  the  expression  for  the  change  of  intercept  of  the  Law  of  the  Wall 
value  of  D  for  sand  roughness 

parameter  in  Han’s  model 

flux  vectors  in  iransfomied  coordinates 
transition  function  from  smooth  to  rough 
height  of  roughness  element 
height  of  equivalent  sand  roughness 
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height  of  sand  roughness 

Clauser’s  constant  =  0.0168 
mixing  length 
free-stream  Mach  number 

streamwise  spacing  of  roughness  elements 
pressure 

vectors  of  dependent  variables 
parameter  in  Han's  model 

Reynolds  number 

reference  area  for  skin  friction,  area  of  smooth  surface  before  adding  on  the 
roughness 

total  frontal  area  of  roughness  elements 
reference  area  for  axial-force  coefficient 

viscous  source  temi 
velocity  in  mainstream  direction 
velocity  at  the  edge  of  the  boundary  layer 
friction  velocity 

longitudinal  distance  from  nose 
distance  from  the  wall 
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r 

K 

X 

A 


P 

CO 


=  ratio  of  specific  iieats 
=  Klebanoff  intermittency  factor 
=  Von  Karman  constant 

=  roughness  densiis  parameter,  defined  in  TABLE  1 
=  composite  roughness  desity  parameter,  defined  in  TABLE  1 
=  turbulent  viscosity 
=  kinematic  viscosity 
=  transfomied  coordinates 
=  density 
=  vorticity 


Superscripts 

i 


o 


r 


s 


inner  layer 
outer  layer 
nnigh 
smooih 

Dalle  Donne  density  parameter 
wall  parameters 
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Subscripts 

a 

f 

h 

mu 

ri 

s 
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density  parameter  base  on  width 

density  parameter  based  on  frontal  area 

density  parameter  based  on  height 

math  wave 

Stan  of  roughness 

new  density  parameter 

wail  conditions 


50 


DISTRIBUTION  LIST 


No. 

Copies  Organization 


No. 

Copies  Organization 


12  Administrator 

Defense  Technical 
Information  Center 
ATTN;  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 

1  HQDA  (SARD-TR) 

Washington,  DC  20310-0001 

1  Commander 

US  Army  Material  Command 
ATTN;  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

2  Commander 
Armament  RD&E  Center 
U.S.  Army  AMCCOM 
ATTN:  SMCAR-MSI 
Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

Armament  RD&E  Center 
U.S.  Army  AMCCOM 
ATTN:  SMCAR-TDC 
Picatinny  Arsenal,  NJ  07806-5000 

1  Commander 

Benet  Weapons  Laboratory 
Armament  Ri:D  Center 
US  Army  AMCCOM 
A'^TN:  SMCAR-LCD-TL 
Watcrvlict.  NY  12189-4030 


1  Director 

US  Army  Aviation  Re.search  and 
Technology  .A.ctivity 
Moffett  Field,  CA  94035-1099 


1  Commander 

US  Army  Communications  - 
Electronics  Command 
ATTN;  AMSEL-ED 
Fort  Monmouth,  NJ  07703-5022 

2  Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-RD 
B.  Walker 

Redstone  Arsenal,  AL  35898-5000 

1  Commander 

US  Army  Missile  Command 
ATTN:  "aMSMI-AS 
Redstone  .\rsenal,  AL  35898-5000 

1  Commander 

US  Armv  Laboratory  Command 
ATTN:  'aMSLC-DL 
Adelphi,  MD  20783-1145 

1  Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  SMCAR-ESP-L 
Rock  Island,  IL  61299-5000 

1  Commander 

US  A  rniy  Aviation  Systems  C<;inrnand 
ATTN:  AMSAV-DACL 
4300  Goodfcllow  Blvd 
St  Louis,  MO  63120-1789 


51 


DISTRIBUTIOxN  LIST 


No. 

Copies  Organization 

1  Commander 

US  Army  Teink  .Automotive 
Command 

ATTN:  AMSTA-TSL 
Warren,  MI  48397-5000 

1  Director 

US  Army  TRADOC  Analysis  Center 

ATTN:  ATAA-SL 

White  Sands  Missile  Range, 

NM  88002-5502 

1  Dir,  Combat  Development 

US  Army  Infeuitry  School 
ATTN:  ATSH-CD 
Fort  Benning,  GA  31905-5GC0 

1  AFWL/SUL 

Kirtland  AFB,  NM  87117-5800 

2  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  R44  (Dr.  F.  Priolo) 

Code  R44  (Dr.  A.  Wardlaw) 
K-24,  Building  402-12 
White  Oak  Laboratory 
Silver  Spring,  MD  20903-5000 


No. 

Copies  Organization 

5  Director 

National  Aeronautics  £uid  Space 
Administration 
Ames  Research  Center 
ATTN:  MS-227-8,  L.  Schiff 
MS-258- 1,  T.  Holst 
MS-258- 1,  J.  Steger 
MS-258-1,  D.  Chaussee 
MS-258- 1,  M.  Rai 
MS-229- 1  M.  Rubesin 
Moffett  Field,  CA  94035 

5  Commander 

U.S.  .Army  .Armament  Research. 
Development  Sz  Engineering  Center 
ATTN:  SMCAR-TSS 

SMCAR-LCA-F 
J.  Grau 

R.  Kline 

S.  Kahn 
H.  Hudgins 

Picatinny  Arsenal,  NJ  07806-5000 

1  .Air  Force  Armament  Laboratory 

ATTN:  AF.ATL/DLODL 
Eglin  AFB,  FL  32542-5000 


3  Director 

National  .Aeronautics  and  Space 
.Administration 
Langley  Research  Center 
ATTN:  Tech  Library 

Mr.  P.  J.  Bobbitt 
Mr.  D.  M.  Eushnell 
Dr.  M.  J.  Hemsch 
Langley  Station 
Hamjrton,  VA  236G5 


United  States  Military  .Academy 
Department  of  Mechanics 
.ATTN:  LTC  .Andrew  L.  Dull 
LTC  Thomas  Kiehne 
West  Point,  NY  1090G 


rj2 


DISTRIBUTION  LIST 


No. 

Copies  Organization 


No. 

Copies  Orsanization 


1  Massachusetts  Institute  of 
Technology 
ATTN:  Tech  Library 
77  Meissachusetts  Avenue 
Cambridge,  MA  02139 

1  Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  Dr.  F.  Moore 
Dahlgren,  VA  22448 

1  Air  Force  Amieiment  Laboratory 

ATTN:  AFATL/FXA 
(Stephen  C.  Korn) 

Eglin  AFB,  FL  32542-5434 

1  AEDC 

Calspan  Field  Semce 

ATTN:  MS  fiOO  (Dr.  John  Benek) 

AAFS,  TN  37389 

1  Virginia  Polytechnic  Institue 

&  State  University 
.\TTN:  Dr.  Clark  H.  Lewis 
Department  of  Aerospace 
Ocean  Engineering 
Blacksburg,  VA  24061 

1  University  of  California, 

Davis 

Department  of  Mechaniccd 
Engineering 

ATTN-  Prof.  H.A.  Dwyer 
Davis,  C.A  95616 

1  Pennsylvania  State  University 

Department  of  Aerospace 
Engineering 

ATTN:  Dr.  G.  S.  Dulikravich 
University  Park,  PA  16802 


1  University  of  Illinois 

at  Urbana  Champaign 
Department  of  Mechaniced  and 
Industrial  Engineering 
Urbeina,  IL  61801 

1  University  of  Mciryland 

Department  of  Aerospace  Engr. 
ATTN:  Dr.  J.  D.  Anderson,  Jr. 
College  Park,  MD  20742 

1  University  of  Notre  Dame 

Department  of  Aeronautical 
eind  Mechanical  Engineering 
ATTN:  Prof.  T.  J.  Mueller 
Notre  Dame,  IN  46556 

1  University  of  Texas 

Department  of  Aerospace 
Engeering  and  Engineering 
Mechanics 

ATTN:  Dr.  D.  S.  Dolling 
Austin,  Texas  78712-1055 

1  University  of  Delaware 

Department  of  Mechanical 
Engineering 

.ATTN:  Dr.  John  Meakin 
Chairman 
Newark,  DE  19716 

1  University  of  Florida 

Department  of  Engineering  Sciences 
College  of  Engineering 
ATTN:  Prof.  C.  C.  Hsu 
Gainesville,  FL  32611 


.53 


DISTRIBUTION  LIST 


No. 

Copies  Organization 

2  USAF  Wright  Aeronautical 

Laboratories 
ATTN:  AFWAL/FIMG 

Mr.  Norman  E.  Scaggs 
Dr.  J.  Shang 
WPAFB,  OH  45433-6553 

1  Grumann  Aerospace  Corporaton 

Aerophysics  Research  Department 
ATTN:  Dr.  R.  E.  Melnik 
Bethpage,  NY  11714 

1  Advanced  Technology  Center 

Arvin/Calspan 
Aerodynamics  Research  Department 
ATTN:  Dr.  M.  S.  Holden 
P.  O.  Box  400 
Buifalo,  NY  14225 

1  Analytical  Methods,  Inc. 

ATTN:  Mr.  F.  A.  Dvorak 
2047  152nd  Ave.,  N.E. 

Redmond,  WA  98052 

2  David  Tavlor  Research  Center 
ATTN:  br.  P.  S.  Granville 

Dr.  de  los  SEmtos 
Bethesda,  MD  20084 

1  Physictd  Science  Inc. 

ATTN:  Dr.  M.  L.  Finson 
Research  Park 
Andover,  MA  01810 

2  Director 

Sandia  National  Laboratories 
ATTN:  Dr.  \V.  Oberkampf 
Dr.  F.  Blottner 
Division  1636 
Albuquerque,  NM  87185 


No. 

Copies  Organization 

1  The  University  of  Arizona 

Aerospace  Engineering  Department 
ATTN:  Prof.  I.  Wygnanski 
Tucson,  AZ  857-1 

1  Applied  Technology  .Associates 

ATTN:  Mr.  R.  J.  CaN^alleri 
P.O.  Box  19434 
Orlando,  FL  32814 

1  United  Technologies  Corporation 

Chemical  Systems  Division 
ATTN:  Mr.  A.  L.  Holzman 
P.O.  Box  50015 
600  Metcalf  Road 
San  Jose,  C.A  95150-0015 

1  Ford  Aerospace  and 

Communications  Corporation 
Aeronutronnics  Division 
ATTN:  Charles  White 
Bud  Blair 
Ford  Road 

Newpoint  Beach,  CA  92658 

1  Honeywell  Inc, 

ATTN:  Wilford  E.  Martwick 
Ken  Sundeen 

600  Second  Street,  North  East 
Hopkins,  MN  55343 

1  University  of  Cincinnati 

Department  of  .Aerospace  Engineering 
ATTN;  Prof.  Stanley  Rubin 
Mail  Location  70 
Cincinnati,  OH  45221 


54 


DISTRIBUTION  LIST 


No. 

No. 

Copies 

Organization 

Copies 

Organization 

1  Illinois  Institute  of  Technology 
College  of  Engineering 
Fluid  Dynamics  Research  Center 
ATTN:  Dr.  Mukund  Acharya 
IIT  Center 
Chicago,  IL  60616 


Aberdeen  Proving  Ground 

Director,  USAMSAA 
ATTN;  AMXSY-D 

AMXSY-MP,  H.  Cohen 

Commander,  USATECOM 
ATTN:  AMSTE-TO-F 

Cdr,  CRDC,  AMCCOM 
ATTN;  SMCCR-MU 

SMCCR-RSP-A 

SxMCCR-SPS-IL 


55 


